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ABSTRACT:  The  formula  of  arrival  angles  of  spatially 
partial  coherent  light  beams  propagating  through 
turbulent  atmosphere  is  derived  from  Markovian  process 
approximation,  which  can  be  used  in  both  weak  and  strong 
intensity  fluctuations.  We  also,  have  derived  the  beam 
wander  angle  and  arrival-angle  fluctuation  variance  of  the 
light  beam  reflected  by  a  mirror  when  the  beam  propagating 
through  the  turbulent  atmosphere.  The  reflecting 
"amplification"  characteristics  of  the  beam  wander  angle 
and  the  arrival  angles  of  light  beams  propagating  along 
a  folded  path  is  discussed,  which  passes  the  distance  twice 
between  the  transmitter-receiver  and  the  mirror  reflector  in 
weak  and  strong  intensity  fluctuations. 

1.  Introduction 
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Due  to  the  needs  of  laser  range  finding  and  atmosphere 
detection  with  laser  radars,  researchers  have  paid  growing 
attention  to  the  statistical  properties  of  various  parameters  of 
beams,  which  twice  pass  through  a  randomly  fluctuating  medium 
with  a  particular  volume.  The  theories [1-8]  and  experiments  [2] 
now  available  suggest  that  in  the  turbulent  atmospheres  with  the 
same  conditions,  the  fluctuation  properties  of  various  parameters 
of  reflected  backward  waves  are  different  from  the  statistical 
properties  of  optical  waves  that  directly  pass  through  the  same 
path . 


In  fact,  the  statistical  properties  of  parameters  of 
backward  waves  are  associated  with  the  geometric  shape  of 
reflectors:  some  reflectors  (such  as  plane  mirror)  can  amplify 
the  deformation  of  optical  waves,  while  others  can  otherwise 
compensate  for  the  deformation  of  optical  waves. 

Consequently,  it  would  be  useful  to  carefully  study  the 
reflection  effects  of  various  reflectors,  and  to  take  advantage 
of  these  effects  in  reducing  or  amplifying  the  atmospheric 
turbulence  effects,  and  to  increase  signal-noise  ratio  in 
measuring  the  atmospheric  parameters. 

However,  up  till  now,  the  discussion  of  reflected  beam 
propagation  direction  fluctuation  has  focused  on  weak  intensity 
fluctuation  alone  in  most  cases[l] .  Although  V.P.  Lukin[2] 
discussed  the  amplification  properties  of  image  trembling  in  the 
receiving  telescope  over  a  longer  propagation  distance,  he  failed 
to  touch  upon  strong  turbulence. 

Through  the  Markov  approximation  and  optical  wave  composite 
amplitude  normal  distribution  approximation,  this  paper  discusses 
the  beam  arrival-angle  fluctuation  while  directly  propagating 
through  the  turbulent  atmosphere,  the  wander  of  reflected 
backward  wave  beam  at  the  receiving  surface  as  well  as  the 
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statistic  properties  of  the  arrival-angle  fluctuation. 

2.  Arrival-Angle  Fluctuation! 

The  experiments  available  indicate  that  the  fluctuations  in 
the  beam  root-mean-square  wander  angle  and  arrival  angle  are  both 
lower  than  10"  [9]  .  Based  on  this,  the  beam  arrival-angle 
fluctuation  problem  can  be  discussed  using  the  paraxial  optics 
approximation.  Under  the  paraxial  optics  approximation,  the 
random  beam  arrival  angle  ag  at  the  receiving  surface  can  be 
derived  as  [8] 

R)t(R)^,  ,  i ) 

where  is  refractivity  fluctuation;  R  is  radius  vector  within 
the  receiving  surface:  t  (R)  =tjexp  (-R‘/a^-  [10]  ;  t^  is  transmittance; 
as  is  the  radius  of  the  receiving  hole;  I=u($,R)u  (£,R)  is  light 
intensity  within  the  plane;  is  horizontal  gradient  operator 

m 

Po-JJ<l(<r,  R)>KR)<PB, 

r  »  .-*•  •  • 

The  beam  is  propagating  along  the  x  axis. 

Assume  that  n,  is  the  random  quantity  of  Gaussian 
distribution  and  that  it  satisfies  the  delta  correlation 
condition,  then  we  can  make  a  systemwide  integration  averaging  of 
the  square  of  Eq.  (1),  and  apply  the  Markov  approximation  [7 , 9] . 

Considering  the  arrival-angle  fluctuation  is  caused  by  the 
optical  wave  phase  fluctuation [2] ,  we  can  simplify  it  using  the 
optical  wave  composite  amplitude  normal  distribution 
approximation [11] .  Then,  through  computation,  we  can  derive  the 
arrival-angle  fluctuation  variance  of  partial  coherent  beams  in 
the  space  of  the  light  source 


3 


<4w 

2(F+ 


(2) 


where 


r  «)  -  [3(i-|)+ (4)’ + (r*+r-)f + (i+-|r>- + r-]r. 

-*\e) -<*£(.£)+**,- 


(3) 


fr/VT:\  f 

The  f  in  Eq.  (2)  is  three-dimensional  turbulent  spectrum 
density;  K  is  number  of  waves  in  space;  ps  is  the  coherence 
length  of  partial  coherent  spherical  waves;  a5  is  space  coherence 
length  of  light  source;  afl  is  the  spatial  coherence  length  of 
light  source;  ac  is  the  equivalent  radius  of  light  source;  k  is 
the  number  of  optical  waves;  F  is  the  curvature  radius  of  beam. 
When  the  light  source  is  a  spherical  wave  or  plane  wave,  Eq.  (2) 
retrogrades  to 

<c£>— d^f 


I  4 ) 


+iexp(-- 


2F 


where  ot^a. .  Eqs .  (2)  and  (4)  are  the  major  results  in  this 


section . 


In  weak  intensity  fluctuation,  when  the  light  source  is  a 
plane  wave  and  Kolmogorov  turbulence  spectrum  is  applied, 
I(C,R)~I(0,R)  [9]  .  Eq.  (4)  presents  <o(‘c>k3<o2:>  =  5 . 76c2nx)  2a,./  ^ , 
which  is  in  agreement  with  the  result  obtained  by  Chiba [4] . 
Here,  <a„>  is  beam  wander  angle  variance. 
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experimental 


Fig.  1.  Comparison  of  the  theory  with  the  experiment 

If  the  Karman  spectrum  or  exponential  spectrum [9]  is  applied 
in  Eqs ,  (2)  and  (4),  they  are  also  suitable  for  strong  intensity 

fluctuation,  which  is  confirmed  by  the  agreement  between  the 
numerical  result*  of  Eq.  (4)  in  Fig.  1  (shown  by  the  imaginary 
line,  and  Karman  spectrum  is  applied  to  <pp  and  the  experiment  [5]  . 
The  solid  line  in  Fig.  1  corresponds  to  the  slight  perturbation 
approximation  theory  [5].  The  result  of 

Fig.  1  also  indicates  that  the  approximation  used  in  deriving 
Eqs.  (2)  and  (4)  is  feasible. 

*  Its  conditions  are:  (x/F)=-l;  Kg=10.4m2;  hg=2m;  v=0.3;  x=1750m; 
o2a=1.06c2nx(2at)  ;  A.=  0.5pm. 


3.  Wander  Fluctuation  of  Reflected  Beam 
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When  a  beam  is  propagating  along  a  folded  optical  path  of 
turbulent  atmosphere  with  a  plane  reflector,  whose  surface  is 
much  larger  than  the  beam  cross  section,  as  the  reflector,  the 
optical  wave  composite  amplitude  u  can  meet  the  following 
parabolic  equation[7]: 

4#»+4WTiC»#  K)v, 

u(p,  *)-*(*), 

where  „  _  (^(x,  R),  x<L,  4 

K)'-\^(2L-x,  K),  L<x<2L, 

An  is  Laplace's  operator;  ufl  (R)  is  optical  wave  composite 
amplitude  at  the  transmission  end. 

The  backward  wave  composite  amplitude  ur  in  L<x£2L  region 
can  satisfy  the  following  equation: 

t  , 

dx  ~  +  X,  R)Ur, 

(6) 

R)-U,(L,  R)a 

where  ur(L,R)  is  the  optical  wave  composite  amplitude  incident  at 
the  reflecting  surface. 


r.Crt 


Fig.  2  Light  beam  propagation  path 

Fig.  2.  Light  beam  propagation  path 

Assume  that  the  center  of  transmitted  beam  is  overlapping 
the  center  of  the  reflected  beam  at  the  receiving  end  with  a 
transmission  optical  path  as  shown  in  Fig.  2.  According  to  the 
beam  wander  theory,  let  £=x,  then  in  the  0<5<L  region,  the  vector 
of  the  beam  "gravity  center"  position  is [8] 
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R)0  '  (7) 

While  in  the  L<§<2L  region,  the  vector  of  the  beam  "gravity 
center"  position  is 

per~'k&  (2i  “  R)V^2L-Z,  J?)0  ( 8  } 

By  ignoring  the  atmospheric  absorption  and  the  backward 
scattering  effect  of  turbulence,  P'0=Pg[2]  can  be  derived  through 
total  beam  flux  conservation.  Then,  by  introducing  this  relation 
in  Eq.  (8),  and  making  variable  replacement  of  £'=£-L, 

Eq.  (8)  will  be  converted  to 

•• 

Pcr-~\liL ~ R)V^L-i,  R)d>R0  ( g } 

The  beam  "gravity  center"  random  wander  angle  at  the 
receiving  end  can  be  calculated  through  superimposition  of  pc/L 
and  Pcr/L  as  follows: 

‘ (L - fidf jj>B(;r(f,  J8)V^«,  J!)  (10) 

R)V^(L-S,  *)].  ..  .. 
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Similar  to  the  derivation  of  Eq.  (2) ,  by  taking  advantage  of 
the  delta  correlation  property  of  nj[7],  the  wander  angle 
fluctuation  variance  of  the  reflected  beam  from  the  plane  ! 

reflector  can  be  derived  j 

mk  m 

X  <!(£,  RJItf,  R») > + <<rfr (L) >+ 2<<t?,  (L) >, 

where  the  first  term  is  the  commonly  used  single  trip  beam  wander 
angle  fluctuation  variance  <oc(L)>;  <ocr(L)  is  the  wander  angle 
fluctuation  variance  of  the  reflected  backward  wave  propagating 
over  a  distance  L;  2<oics(L)>  is  the  contribution  of  the  increased 
phase  fluctuation  caused  by  the  interference  between  the  incident 
light  and  reflected  light  to  the  beam  wander.  Their  specific 
expressions,  respectively,  are: 

( 12 ) 

x<I,(l-£,  RJI, (l-£,  R,)>expiiK-(R1-Ra)l, 

(13) 

x <I(£  RJIrtf,  «a)>exp[iK.  («!-«,)] e 

Obviously,  <o%A  ( 2L)  >><o2c  ( 2L)  >  .  Here,  <o2c(2L)>  is  the 
wander  angle  fluctuation  variance  of  the  beam  propagating  through 
the  2L  turbulence  layer. 


In  weak  Kolmogorov  turbulence  fluctuation  as  well  as  in 
collimated  and  focused  beam  propagation,  the 

I  ( 5  ,  R)  ~I  ( 0  ,  R)  »Ir  ( L-£  ,  R )  approximation  can  be  made[9].  Under  this 
approximation ,  .<^(i)>-2<£r*.(L)>- 

<cr2(£)>  -  [0 . 099*?  r  (1/6)  c2/2*/2"  ]  L  (2ao)  -1/*  -  (1/2)  <<r,  (2L)  >  can  be  derived 

from  Eqs .  (12)  and  (13).  By  inserting  this  relation  in  Eq.  (11), 

the  following  can  be  obtained: 


<<r**(2L)>-A<<r#(2L)>. 


(14) 
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This  result  conforms  to  the  following  result:  a  finite  plane 
wave,  derived  through  directly  and  approximately  applying  the 
logarithmic  normal  distribution  of  optical  wave  composite 
amplitude  in  the  wander  fluctuation  variance  and  defining 
<P=’>-  Ru  Rt,  Hi,  RJ/Pl  ,  is  reflected  from 

the  plane  reflector  and  then  wander  amplified  by  3/2  times [1]. 
However,  our  conclusion  is  also  applicable  to  the  focused  beam 
propagation. 

In  strong  intensity  fluctuation,  we  now  discuss  the  wander 
amplification  problem  based  on  the  analysis  of  the  optical  wave 
mutual  coherence  function  Ar( P,  R)  presented  in 

ref erence [8] .  According  to  the  commutation  principle,  when  a 
beam  propagating  through  a  turbulence  layer  eC^'  is 

reflected  from  a  plane  reflector  located  at  x=L,  the  mutual 
coherence  function  of  the  reflected  optical  wave  is [8] 

xo-(o;  «„-!,*  i,  s+I*  l,  R’+~p') 

***(*'.*—§*  L-  (i5) 

where  u(p)  is  the  optical  wave  composite  amplitude  incident  at 
the  turbulence  layer  x=0;  u2L(p)  is  the  reflected  wave  composite 
amplitude  at  x=2L;  G(0,pa;  x,p)  is  the  spherical  wave  composite 
amplitude  in  the  turbulence  medium,  which  can  satisfy  the 
following  equation: 

2ik-^-+AMG-‘2t>n1(aJ  p)G,  G(0,  pw  m,  p)-&(p*-p)0 

Suppose  the  turbulence  fluctuation  is  sufficiently  strong, 
and  after  the  spherical  wave  passes  through  the  turbulence  layer 
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with  a  distance  L,  its  composite  amplitude  can  meet  the  normal 
distribution  statistics,  then  the  following  can  be  derived  from 
the  statistic  property  of  the  norma!l  distribution: 


ru(p,  R) -r*(L,  p,  R)+r»(L,  P,  R), 

X (*-*.) -D.O',  P)-D.(fi',  V-C)], 

w.  »■ ») -(is- (*.+p-i-  ?>•[*•  -(pr-^  />)] 
—D»{p',R—R<>+pr)t 


(16) 


where  Ds  is  phase  structure  function[8];  r  jr  is  the  mutual 
coherence  function  of  the  optical  wave  that  directly  passes  a 
distance  2L;  r2r  is  the  long-distance  correlation  function  of  the 
reflected  field,  which  is  related  to  the  pre-reflection  and  post¬ 
reflection  beam  propagation  through  the  same  inhomogeneous 
medium.  When  the  turbulence  fluctuation  is  sufficiently  strong 
to  make  the  beam  turbulence  spread  much  larger  than  diffraction 
spread,  <I2r  (LfR)  )  >»<Ilr(L,R)  >,  i.e.  <I2L  (L,R)  >~<Ilr  (L,R)  >  [8]  can  be 
derived  from  Eq.  (16).  Here,  <Ijr  2r  (L ,  R)  >  =  r2r  2r  (L ,  0  ,  R)  .  . 

The  foregoing  discussion  indicates  that  when  turbulence 
spread  is  much  larger  than  diffraction  spread,  the  reflection 
effect  of  the  plane  reflector  r 2r  is  fairly  weak,  and 
accordingly,  the  contribution  of  2<o^cs(L)>  to  <o^(2L)>  is  fairly 
small,  which  will  be  further  explained  in  the  following 
discussion.  In  the  strong  intensity  fluctuation,  by  using 
optical  wave  composite  amplitude  normal  distribution 
approximation,  the  light  intensity  mutual  correlation  function  of 
the  reflected  field  is  Ri))>a 

Ri)XI*l(£,  R»)>+\rM  Ru  R*)\a#<i*.(.t,  RJXiM  «.)> 

[9]  .  Since  it  is  known  from,  the  foregoing  discussion  that 
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!<**(&  **)>»<!*(£  *l)Ilr(£,  *»)>  ,  then  from  the  wander 

formula [7] : 

<a*A(2i)>-4 »  (W)*#}/  iPKK'MX)  it<PR, 

x  exp[iK‘ (Ri— Ri)IiL(.£,  R»)y, 

it  is  known  that  <o^ci  ( 2L)  >  =  <o^c  ( 2L)  >  .  This  suggests  that  when  the 
ratio  between  turbulence  spread  and  diffraction  spread  increases 
with  the  increase  in  the  optical  wave  propagation  distance  or  in 
turbulence  intensity,  the  wander  amplification  effect  of  the 
plane  reflector  decreases. 


4.  Arrival-Angle  Fluctuation  of  Reflected  Beam 

Based  on  the  above  discussion,  at  the  receiving  end,  the 
random  arrival  angle  of  the  beam  wavefront  reflected  from  the 
plane  reflector  placed  in  the  middle  of  the  optical  path  can  be 
expressed  as 

dEoL-ipi+ien,  d7) 

da  dx  dx 

where  (dpcl/dx)  is  the  wavefront  fluctuation  caused  by  turbulence 
in  the  0^|;<L  region;  (dp^/dx)  is  the  wavefront  fluctuation 
caused  by  turbulence  in  the  L<$<2L  region. 

Similar  to  the  derivation  of  Eq.  (2),  the  reflected  beam 
arrival-angle  fluctuation  variance  can  be  derived  from  Eqs .  (1) 

and  (17)  as  follows: 

'  + *(■&■■&)■ 

where  <o^c(2L)>  is  the  arrival-angle  fluctuation  variance  of  the 
beam  directly  propagating  through  the  turbulence  layer  2L;  the 
second  term  of  Eq.  (18)  is  the  contribution  of  optical  wave 
interference  before  and  after  reflection  to  the  arrival-angle 
fluctuation,  its  specified  expression  being 
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(19) 


'  '  '  ■  x<I((,  *)I,(f,  *,)><(*.)<(*.). 

Ill  weak  intensity  fluctuation,  when  a  collimated  and  focused 
beam  serves  as  the  light  source,  the  I(£,  K)&I,(£, 
approximation  can  be  made.  Under  this  approximation, 

.  From  this  relation  and 
Eq.  (17) ,  the  following  can  be  obtained: 

'  <«?,(2i)>-2<a?(2L)>0  .  (20) 

This  result  conforms  to  the  result  acquired  using  other 
techniques  as  shown  in  ref erences  [1]  and  [10].  Similar  to  the 
discussion  in  the  previous  section,  when  the  turbulence  spread  is 
much  larger  than  diffraction  spread,  <orcr(2L)>  becomes 
<a^c(2L)>.  This  shows  that  not  only  the  wander  amplification 
effect  of  the  plane  reflector  decreases  with  the  increase  in 
propagation  distance,  but  also  the  reflected  beam  arrival-angle 
fluctuation  amplification  also  decreases  when  the  turbulence 
intensity  C^n  increases  so  that  the  turbulence  spread  is  much 
larger  than  the  diffraction  spread. 

If  a  reflection  kernel  function  is  introduced  in  the 
foregoing  discussion [1] ,  then  it  can  be  known  through  similar 
discussion  that  a  triangle  prism  can  reduce  the  reflected  beam 
wander  angle  and  arrival-angle  fluctuation  in  the  weak  intensity 
fluctuation . 

5.  Conclusions 

With  the  "Markov  approximation",  this  paper  discussed  the 
arrival-angle  fluctuation  of  partial  coherent  beams  in  the 
turbulent  atmosphere,  and  derived  simple  expressions  applicable 
for  the  entire  fluctuation.  Also,  we  derived  formulas  of  beam 
wander  angle  and  arrival-angle  fluctuation  at  a  folded  path, 
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which  prove  to  be  even  more  universal  compared  with  the  theories 
available  now.  Finally,  we  discussed  the  reflection 
amplification  property  of  the  beam  propagation  direction 
fluctuation  in  the  strong  intensity  fluctuation,  and  derived  the 
amplification  rate  of  the  wander  angle  and  arrival  angle 
fluctuation  of  the  collimated  and  focused  beam  reflected  from  the 
plane  reflector  in  the  weak  intensity  fluctuation. 
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